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Abstract. Hadron spectra from the first year RHIC run are shown to be excellently 
reproduced by hydrodynamic calculations. We argue that in particular the elliptic flow 
data provide strong evidence for early thermalization at RHIC, at energy densities well 
above deconfinement, but that the phenomenologically extracted short thermalization 
time scale of less than 1 fm/c provides a serious challenge for theory. The HBT radii 
from the hydrodynamic calculations agree only qualitatively with the data, showing 
significant quantitative discrepancies. It is argued that this points to a still incomplete 
understanding of the freeze-out process at RHIC. 

Keywords: Relativistic heavy-ion collisions, hydrodynamics, thermalization, two- 
particle correlations 

PACS: 25.75. -q, 25.75.Ld, 25.75.Gz, 24.10.Nz 



1. Elliptic flow as an early QGP signature 

The quark-gluon plasma (QGP) is a thermalized system and, as such, has thermal pressure. 
If the QGP is created in a heavy-ion collision, this pressure acts against the surrounding 
vacuum and causes a rapid collective expansion ("flow") of the reaction zone, the "Little 
Bang". In non-central collisions the initial overlap region of the two nuclei is elliptically 
deformed in the transverse plane, resulting in anisotropic pressure gradients. These cause 
a more rapid expansion into the reaction plane than perpendicular to it, resulting in an 
anisotropy of the final transverse momentum distribution called elliptic flow [|l|]. Elliptic 
flow is quantified by the second harmonic coefficient y2{y,P±',b) of a Fourier expansion in 
of the measured hadron spectrum dN/ {dy p±dp± d(^p) [ ^. At midrapidity in collisions 
between equal-mass nuclei it is the lowest-order non-zero Fourier coefficient. 

Since individual nucleon-nucleon collisions produce azimuthally symmetric spectra, 
any final state momentum anisotropics must be generated dynamically during the nuclear 
reaction. A non-zero value for the elliptic flow is thus an unambiguous sign for reinterac- 
tions among the produced particles, and its magnitude characterizes their strength. Micro- 
scopic transport calculations [ |[ ^ show a monotonic dependence of V2 on the opacity 
(density times scattering cross section) of the produced matter which is inversely related 
to its thermalization time. Elliptic flow also requires the existence of an initial spatial 
anisotropy of the reaction zone, by colliding either deformed nuclei at ^7 = or spherical 
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nuclei at b^O. The transport calculations strongly suggest that, for a given initial spatial 
anisotropy e^, the maximum momentum-space response V2 is obtained in the hydrodynamic 
limit which assumes a thermalization time which is much shorter than any macroscopic 
time scale in the system. Any significant delay of thermalization allows the initial spa- 
tial deformation to decay without concurrent build-up of momentum anisotropics, thereby 
reducing the finally observed elliptic flow signal [ |]]. 

This specific sensitivity of the elliptic flow to rescattering and pressure build-up in the 
early collision stages [ ||, ^ (before the spatial deformation and the resulting anisotropics 
of the pressure gradients have disappeared [ ^]) puts V2 on the list of "early signatures" of 
the collision dynamics. In contrast to other early probes (which use rare signals such as 
hard photons and dileptons, heavy quarkonia and jets), V2 can be extracted from the bulk 
of the measured hadrons which are very abundant and thus easily accessible. In fact, the 
elliptic flow measurement in AuH- Au collisions at y/s = 1 30 A GeV [ |8| became the second 
publication of RHIC data and appeared within days of the end of the first RHIC run. 

We present hydrodynamic results for hadronic spectra and elliptic flow at RHIC en- 
ergies. We show that the hydrodynamic approach provides an excellent quantitative de- 
scription of the bulk of the data and fails only for very peripheral AuH-Au collisions and/or 
at high px>l -5— 2GeV/c. That the hydrodynamic approach fails for small initial overlap 
regions or large hadron transverse momenta is not unexpected. What is surprising is that 
the hydrodynamic approach does work for such a wide range of impact parameters and 

™d quantitatively reproduces the momenta of more than 99% of the particles: below 



px = l-5GeV/c the elliptic flow data [ |^ |9|, [10| actually exhaust the hydrodynamically 
predicted [ I [n], |l2[ |l^] upper limit. The significance of this agreement can hardly be 
overstressed, and it poses significant challenges for microscopic descriptions of the early 
collision dynamics. How the system is able to thermalize so fast is the first RHIC puzzle 
which we point out. 

The second puzzle arises from two-particle correlation measurements (HBT interfer- 
ometry) [ [l^, |l^ which constrain the freeze-out distribution in space-time. The same hy- 
drodynamic model which provides an almost perfect description of the momentum-space 
structure of the emitting source (spectra an elliptic flow) predicts a space-time distribution 
at freeze-out which, when translated into HBT radii, does not agree very well with the 
data. We argue in Sec. || that this points to a problem with our understanding of the late 
freeze-out stage of the collision. 



2. Radial and elliptic flow from hydrodynamics 

The natural language for describing collective flow phenomena is hydrodynamics. In the 
ideal fluid (non-viscous) limit used by us, this approach assumes that the microscopic mo- 
mentum distribution is thermal at every point in space and time. This does not require 
chemical equilibrium - chemically non-equilibrated situations can be treated by solving 
separate and coupled conservation equations for the particle currents of individual particle 
species. The assumption of local thermal equilibrium is an external input, and hydrody- 
namics offers no insights about the equilibration mechanisms. It is clearly invalid during 
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the initial particle production and early recattering stage, and it again breaks down towards 
the end when the matter has become so dilute that rescattering ceases and the hadrons 
"freeze out". The hydrodynamic approach thus requires a set of initial conditions for the 
hydrodynamic variables at the earliest time at which the assumption of local thermal equi- 
librium is applicable, and a "freeze-out prescription " at the end. For the latter we use the 
Cooper-Frye algorithm [ |l7|] which implements an idealized sudden transition from per- 
fect local thermal equilibrium to free-streaming. A better algorithm [ [l^, |l2|] switches 
from a hydrodynamic description to a microscopic hadron cascade at or shortly after the 
quark-hadron transition, before the matter becomes too dilute, and lets the cascade handle 
the freeze-out kinetics. The resulting flow patterns [ |l^ from such an improved freeze-out 
algorithm don't differ much from our simpler Cooper-Frye based approach. 

The main advantage of the microscopic freeze-out algorithm [18, |l^ is that it also 



correctly reproduces the final chemical composition of the fireball, since the particle abun- 
dances already freeze out at hadronization, due to a lack of particle-number changing in- 
elastic rescattering processes in the hadronic phase [|l|]. Our version of the hydrodynamic 
approach uses an equation of state which assumes local chemical equilibrium all the way 
down to kinetic freeze-out at rf«125MeV and thus is unable to reproduce the correct 
hadron yield ratios. We therefore adjust the normalization of the momentum spectra for the 
rarer particle species (kaons, protons, antiprotons) in central collisions by hand to repro- 
duce the chemical equilibrium ratios at a chemical freeze-out temperature T^chem ~ 

165 MeV. 

The absolute normalization of the pion spectra is adjusted through the initial energy density 
in central collisions; for non-central collisions no new parameters enter since the centrality 
dependence of the initial conditions is completely controlled by the collision geometry. 

To simplify the numerical task of solving the hydrodynamic equations we analytically 
impose boost invariant longitudinal expansion [ ^ |l]|. This doesn't give up any essential 
physics as long as we focus on the transverse expansion dynamics near midrapidity (the 
region which most RHIC experiments cover best). The hydrodynamic expansion starts at 
time Teq which we fixed by a fit to hadron spectra from the SPS and then extrapolated to 
RHIC initial conditions (for details see [ ^ |l3[|). For each impact parameter the initial 
energy density profile in the transverse plane is calculated from a Glauber parametrization 
using realistic nuclear thickness functions [ |], ^] . The measured centrality dependence 
of the charged particle rapidity density at midrapidity, (dNch/dy) (y = Q), selects certain al- 
lowed combinations of "hard" and "soft" mechanisms for the initial particle production [ 
pT| ]. We here present results for initial conditions at Teq=0.6fm/c calculated from a mix- 
ture of 25% "hard" (binary collision) and 75% "soft" (wounded nucleon) contributions 
[ ^1]] to the initial entropy density (or parton density), with a maximal entropy density 
■Smax ~ 85 /fm^ at the fireball center in central collisions (corresponding to a maximal en- 
ergy density Cmax = 2 1 .4 Ge V/fm"^ and a maximal te mpe rature Ty^ax ~ 

328 MeV). At the 

standard time T = 1 fm/c used in Bjorken's formula [ for estimating the energy den- 
sity from the measured multiplicity density, this corresponds to an average energy den- 
sity (e)(lfm/c)=5.4GeV/fm^ which is about 70% higher than the value reported from 
158AGeV Pb-nPb colHsions at the SPS. (Note that (e) at Teq =0.6 fm/c is nearly twice as 
large!) The corresponding profiles for peripheral collisions are then given by the Glauber 
model [ 0, |T]|. Kinetic freeze-out was forced at 7f = 130 MeV, independent of centrality. 



4 



U. Heinz and P.F. Kolb 



10' 



Q. 

■D 



10" 





It* PHENIX 

* p PHENIX 

• p STAR 




— i:"^ hydro 

- - - p hydro 


most central 





PHENIX 



centrality- 



p^(GeV) 




p^(GeV) 



PHENIX 




> 

a 

H 
Q. 
■D 

H 
Q. 

>s 

z 
■o 

CM 



10 



10" 



10" 



10" 



- 5 % 


. PHENIX 
— hydro 


-5-16% 












centrality i 









Fig. 1. 



p^(GeV) 



0.5 



1.5 



p^(GeV) 



Charged pion, antiproton and positive kaon spectra from central (upper left panel) and 
semi-central to peripheral (other three panels) Au+Au collisions at ^/s= 130 A GeV. The data were 
taken by the PHENIX [ |^ and STAR [ collaborations (the STAR data have slightly different 
centralities than the indicated values from PHENIX). The curves show hydrodynamical calculations 
(see text). 

Fig. [l| shows the (absolutely normalized) single particle ;?x-spectra for charged pions, 
antiprotons and kaons measured in Au+Au collisions at RHIC together with the hydro- 
dynamical results. The latter were normalized in central collisions as described above, 
but their centrality dependence and shapes are then completely fixed by the model. The 
agreement with the data is impressive; for antiprotons the data go out to p±<3 GeV/c, 
and the hydrodynamic model still works within errors! Only for very peripheral collisions 
(impact parameter b> 10 fm) the data show a significant excess of high-px particles at 
p±>1.5 GeV/c. Teaney et al. [ |l^ showed that this excellent agreement requires a phase 
transition (soft region) in the equation of state; without the transition, the agreement is lost, 
especially when the constraints from SPS data and from the elliptic flow measurements 
below are taken into account. 

In the hydrodynamic picture the fact that antiprotons become more abundant than pi- 
ons for p± > 2GeV/c (upper left panel of Fig. |l|) is not surprising at all: it is a simple 
consequence of the strong radial flow at RHIC. For a hydrodynamically expanding thermal- 
ized fireball, at relativistic transverse momenta p± ^ mo all hadron spectra have the same 
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shape [ and at fixed mj^^niQ their relative normalization is given by [giki) / [gjkj) 
(where gij is the spin-isospin degeneracy factor and A-;,; = e^'-'l^ is the fugacity of hadron 
species At RHIC the baryon chemical potential at chemical freeze-out is small, 

(UB/7chem~0.26 [ |7]|, and /Jtc^O; the pl%^ ratio at fixed and sufficiently large mx is thus 
predicted to be larger than 1: {p/%^),n^ =2exp[— (;UB+;U,i)/rchem] ~ 1-5 (where the factor 
2 arises from the spin degeneracy of the p). 




Fig. 2. The elliptic flow coefficient V2(px) for all charged particles (left) and for identified pions 
and protons (right) from 130AGeV minimum bias Au+Au collisions [ ^, |^ ^ ^9||. The curves 
are hydrodynamic calculations corresponding to equations of state with (Q) and without (H) a phase 
transition and (in the right panel) three different freeze-out temperatures (Tf = 1 28 MeV (dash-dotted), 
BOMeV (solid) and 134 MeV (dashed)). 

Figure ^compares the differential elliptic flow V2(px) from minimum bias AuH-Au 
collisions at RHIC [ ^ ^ 29 1 with hydrodynamic calculations. For transverse momenta 
px < 2 GeV/c the data are seen to exhaust the upper limit for V2 obtained from the hydro- 
dynamic calculations. Only for hadrons with px > 2GeV/c (i.e. fewer than 1% of all) vi 
stays below the hydrodynamic upper limit, indicating incomplete thermalization of high- 
px particles. The data also show the hydrodynamically predicted mass-dependence of V2 [ 
p4[ ] (right panel). Similar conclusions can be drawn from a plot of the px -integrated elliptic 
flow vs. centrality [ 13, only for large impact parameters b>lfm the measured V2 



remains significantly below the hydrodynamic prediction, presumably indicating a lack of 
early thermalization when the initial overlap region becomes too small. 

The excellent agreement with hydrodynamics becomes even more impressive after you 
begin to realize how easily it is destroyed: As stressed in Sec. |l], it requires the build-up of 
momentum anisotropies during the very early collision stages when the spatial anisotropy 
of the reaction zone is still appreciable, causing significant anisotropies of the pressure 
gradients. A delay in thermalization by more than about 1 fm/c (2 fm/c) dilutes the spatial 
anisotropy and the hydrodynamically predicted elUptic flow coefficient by 10% (25%) [ |[ 
which is more than is allowed by the data. Parton cascade simulations with standard HI- 
JING input generate almost no elliptic flow and require an artificial increase of the opacity 
of the partonic matter by a factor 80 to reproduce the RHIC data [ ^. Hadronic cascades 
of the RQMD and URQMD type (in which the high-density initial state is parametrized 
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by non-interacting, pressureless QCD strings) predict [ ^ too little elliptic flow and a 
decrease of V2 from SPS to RHIC, contrary to the data. One can get close to the data by 
forcing the strings to "melt" into partons which then rescatter with large cross sections [ 
^H ]. What causes these large cross sections is still a mystery. 

The elliptic flow is self-quenching [ ||]: it makes the reaction zone grow faster along 
its initially short direction and thus eventually eliminates its own cause. As the spatial 
deformation of the fireball goes to zero, the elliptic flow saturates [ ||]. The saturation 
time scale times c is of the order of the transverse size of the initial overlap region (at 
lower energies it is a bit longer, see Figs. 7,9 in [ |]). At RHIC energies and above, the 
time it takes the colhsion zone to dilute from the high initial energy density to the critical 
value for hadronization is equal to or longer than this saturation time: most or all of the 
elliptic flow is generated before any hadrons even appear! It thus seems that the only 
possible conclusion from the successful hydrodynamic description of the observed radial 
and elliptic flow patterns is that the thermal pressure driving the elliptic flow is partonic 
pressure, and that the early stage of the collision must have been a thermahzed quark-gluon 
plasma. 



3. The RHIC HBT puzzle 

Hydrodynamics not only predicts the momenta of the emitted hadrons, but also the spa- 
tial structure of the hadron emitting source at freeze-out. Bose-Einstein (a.k.a. Han- 
bury Brown-Twiss (HBT)) two-particle intensity interferometry allows to access the rm.s. 
widths of the space-time distribution of hadrons with a given momentum p [ One of 
the interesting questions one can try to address with this tool is whether at RHIC the reac- 
tion zone flips the sign of its spatial deformation between initial impact and final freeze-out, 
as expected at very high collision energy [ ||]. The answer to this question turns out to be 
non-trivial, on two different levels: first, hydrodynamics, at least with the presently im- 
plemented initial conditions and freeze-out algorithm, fails to reproduce even for central 
AuH-Au collisions the measured HBT radii extracted from two-pion correlations [ |l5[ [T^ ]. 
We'll show how and explain why. Second, for expanding systems the HBT radii don't 
measure the entire freeze-out region, but only the effective emission regions ("regions of 
homogeneity") for particles of given momentum [ |3^. For non-central collisions, due to 
the anisotropic transverse flow these can, at least in principle, have a different spatial defor- 
mation than the entire (momentum-integrated) freeze-out region, giving rise to a different 
behaviour of the HBT radii observed at different angles relative to the reaction plane than 



perhaps naively expected [ 33|. 

Fig. ^ shows the HBT radii /?side^ ^out, and R\oaa for central AuH-Au collisions at 
y/s= I30AGeV [ |l5|, 0] compared with hydrodynamic results. The theoretical results 
only account for directly emitted pions; resonance decay pions would add an exponential 
tail to the emission time distribution and could cause a slight increase in the emission du- 
ration, without affecting much the spatial width of the emission region [ ^5[|. The solid 
lines correspond to initial conditions and freeze-out parameters as used for the successful 
description of the momentum distributions in the previous section; clearly, RsHt comes out 
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Fig. 3. HBT radii from a hydrodynamic source compared to RHIC data [ |l5|, The solid lines 
show hydrodynamic results with standard initialization and freeze-out (see text). The dotted lines 
assume freeze-out directly after hadronization at ejec =«crit- The other lines correspond to modified 
initial conditions as described in the text. 

too small whereas /?out and /?iong are both too large in the model. The problem of Rf^i^^ from 
hydrodynamics being too small is well-known from the SPS [ |3^; presumably it is mostly 
due to the sharp Cooper-Frye freeze-out and seems to be at least partially resolved if the 
freeze-out kinetics is handled microscopically within a hadronic cascade [ The latter 
gives a more "fuzzy" spatial freeze-out distribution with larger rm.s. width in the sideward 
direction; it is not so clear that it also resolves the problem with the /Tx -dependence of R^ide 
which is much stronger in the data than in the model. 

On the other hand, "fuzzy" freeze-out [ ^] and the inclusion of resonance decay con- 
tributions [ ^5|] only exacerbate the problems with R]ong and Rout- In high energy heavy-ion 
collisions the longitudinal HBT radius is controlled by the expansion dynamics via the 
longitudinal velocity gradient at freeze-out [ For a boost invariant longitudinal flow 
profile this gradient decreases with time as 1 /x, leading to rather weak gradients (and corre- 
spondingly large values for R]ong) at the typical hydrodynamic freeze-out time of 15 fm/c. 
To reduce ^long, either freeze-out must happen significantly earlier, or the expansion of 
the source must strongly violate longitudinal boost-invariance even close to mid-rapidity, 
somehow allowing for significantly stronger longitudinal flow velocity gradients at late 
times than predicted by the Bjorken profile vi=z/t. 

Noting that decoupling is driven by the strong radial flow developing during the ex- 
pansion stage [ we have tried to force the system to decouple earlier by initiating the 
transverse flow even before Xeq. Even if the system is not yet locally thermalized, it will 
start to develop some transverse collective dynamics, albeit perhaps not as quickly as in 
a hydrodynamic approach. We have tested two extreme assumptions about the transverse 
expansion prior to thermalization [ ^]: in one simulation, shown as the dot-dashed curve 
in Fig. ^ we started the hydrodynamic evolution directly at the parton formation time (for 
which we took the somewhat arbitrary value Tfonn =0.2 fm/c). In another limit (dashed 
lines in Fig. |]), we allowed the partons to stream freely from time Tform to Teq and matched 
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at Teq the first row of the energy momentum tensor to an ideal fluid form, thereby extracting 
an initial transverse flow profile at Xeq. In both cases the resulting transverse flow "seed" at 
Teq = 0.6 fm/c caused the system to expand more rapidly and farther out into the transverse 
direction, freezing out 10-20% earlier. Fig. ^ shows that this helps with both Riong and Rout, 
but not as much as required by the data. 

To get close to the data, we would need to postulate freeze-out directly at the hadroni- 
zation point (dotted lines in Fig. |3]). This radical and entirely unmotivated assumption 
almost removes the discrepancy with R\ong and with the average magnitude of Rom, but it 
is still unable to reproduce the measured strong ATi-depenence of Rout and i^side^ and it stiU 
yields much too small values for Rf^n^ at low K±. Correspondingly, the ratio ^out/^side 
(which is w 1 in the data) is still significantly overpredicted by the model (although much 
less so than by earlier hydrodynamic predictions [ when constraints from RHIC mo- 
mentum spectra were not yet available). 

The problem with the small R^ide cannot be fixed by pre-equilibrium transverse flow, 
either: even though the system then expands to larger values of r, R^^i^e remains essentially 
unchanged (dashed and dash-dotted lines in Fig. H). The homogeneity region simply moves 
farther out without increasing in size. We don't see a way to move closer to the data by 
further modifying the initial conditions, and even the alterations we made to obtain Fig. g 
may turn out to be excluded by the singles spectra and the elliptic flow data (which we 
haven't tested yet). We therefore believe that a resolution to the HBT puzzle must lie in the 
handling of the freeze-out process (although we don't yet know how). 
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Fig. 4. The AT^ -integrated emission function from hydrodynamics, for directly emitted pions with 
rapidity Y = from semicentral Au+Au coUisions (fo = 7fm). The emission function is integrated 
over time and longitudinal coordinate z- The left panel shows contours of constant particle density 
in the transverse plane, the right panel presents cuts through this distribution along the x and y axes, 
showing that particle emission is strongly concentrated at the surface, especially for emission into the 
reaction plane (y = 0). 

Let us further analyze the relationship between Rout and /?side- The formula [ B2p 



^out ^side ~ 



out-.^out)^ - fede--Xside)^) + ^±{it-tf) - 2^j_{(Xo 



t)(f-f)> (1) 
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shows that their difference is controlled by three contributions which can partially com- 
pensate each other There is always a positive contribution from the emission duration 
(5f )" = ((f — f )^) ; for systems with longitudinal boost-invariance freeze-out happens at con- 
stant longitudinal proper time T = Tf, and late freeze-out, which leads to a large longi- 
tudinal homogeneity region 7?iong, thus implies contributions from a wide region along 
this Z — Zf hyperbola, leading to a large variance 5f. To reduce 5f we should therefore 
again freeze out earlier. In hydrodynamics with Cooper-Frye freeze-out [ [l^], the term 
—2^±{{xout—Xout){t—t)) is also positive, since freeze-out happens from the outside inward, 
causing a negative Xom—t correlation. These two positive contributions to ^out^^side 
partially cancelled by a negative geometric contribution (SjCout)^ — (S^side)^: according to 
Figures ^ and ^ hydrodynamics predicts that pion emission is strongly surface dominated, 
especially for pions with non-zero transverse momentum K±, whose emission regions are 
tightly squeezed towards the edge of the fireball where the radial flow is strongest. This 
results in an "opaque source", characterized by a smaller outward than sideward variance 
(5jCout)^ < (5-fside)^ [ Ell- The opacity of the hydrodynamic source is stronger at RHIC 
than SPS energies, but not strong enough to compensate for the two positive contributions 
in Eq. (|l]); this is the reason why we fail to reproduce the measured relation Rout ~ ^side- 
Of course, things would be easier if the Xom—t correlation were positive (i.e. pions at larger 
radial distances were emitted later), but we don't see a good reason why this should be so. 

Let us close with a quick preview of results for the HBT radii for non-central collisions, 
in particular their dependence on the angle <I> of the transverse pair momentum K± relative 
to the reaction plane. Fig. ^ shows our hydrodynamic results for /^jye' ^ouf ^os' ^^'^ ^k)ng 
for b — 7fm AuH-Au collisions at a/s= DOAGeV, plotted as functions of <I> for a number 
of different values of K± ~ \K±\. (Note that "out" and "side" denote the directions parallel 
and perpendicular to in the transverse plane [ ^2[|.) While ^fo^g is almost independent 
of <I>, the three other radius parameters show marked azimuthal dependences of the generic 
form (with all coefficients being positive) 

4de(1>) -<0+'Ri2COs(2<I>), =<0-<2COS(2<I>), 

^os(^)=<2sin(2<I>). (2) 

Although the magnitudes of the coefficients in Fig. ^are quite different (and presumably 
not too trustworthy, given the disagreement with the data for central collisions in Fig. H), 
it is surprising that the signs and phases of the oscillations are identical to those calculated 
and measured at the AGS [ ! At the AGS radial flow effects are thought to be sufficiently 
weak that the oscillations can be interpreted purely geometrically [ reflecting a spatially 
deformed source which is elongated perpendicular to the reaction plane (as is the case for 
the initial overlap region). Fig. ^ shows that at RHIC the reason for these oscillations 
is more subtle: except for small K±, where pions are emitted from the entire fireball, the 
strong radial flow squeezes the emission regions towards the edge of the fireball so that they 
no longer reflect directly the spatial deformation of the momentum-integrated emission 
function shown in Fig. Q However, even though the measured HBT radii have little to do 
with the total widths of the fireball in x and y directions, the out-of -plane deformation of 
the latter imprints itself on the shape of the effective emission regions through the "opacity 
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Fig. 5. Time- and z-integrated emission function for b = 7fm Au+Au collisions at RHIC for 
directly emitted pions with rapidity Y = and fixed transverse momentum K^. Shown are contours 
of constant particle density at freeze-out in the transverse plane, for A'x = in the upper left panel, 
and for A'^ =0.5 GeV and three azimuthal emission angles (0°, 45°, and 90° relative to the reaction 
plane) in the remaining three panels. 

effect", i.e. the fact that the emission is strongly surface-peaked, leading to azimuthal 
oscillations which are in phase with those expected from naive geometrical considerations 
of the entire source. 

Even though the absolute magnitudes of the radii shown in Fig. ^ are affected by 
similar problems as those in Fig. ^ the phases of the azimuthal oscillations agree with 
preliminary STAR data [ and also their amplitudes are at least in the right ballpark. This 
gives rise to some optimism that the HBT puzzle is limited in scope and can be resolved by 
evolutionary rather than revolutionary methods. 
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